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ABSTRACT 

We present new optical spectroscopy for 342 R < 18 galaxies in the Shapley Su- 
percluster (and 198 supplementary galaxies), obtained from 8 hour integrations with 
the AAOmega facility at the Anglo- Australian Telescope. We describe the observations 
and measurements of central velocity dispersion a, emission line equivalent widths and 
absorption line indices. The distinguishing characteristic of the survey is its coverage 
of a very wide baseline in velocity dispersion (90 per cent range a = 40 — 230 km s -1 ), 
while achieving high signal-to- noise ratio throughout (median 60 per A at 5000 A). 
The data quality will enable estimates of Balmer-line ages to better than 20 per 
cent precision even for the faintest galaxies in the sample. Significant emission at Ha 
was detected in ^30 per cent of the supercluster galaxies, including ~20 per cent 
of red-sequence members. Using line-ratio diagnostics, we find that the emission is 
LINER-like at high luminosity, but driven by star-formation in low-luminosity galax- 
ies. To characterise the absorption lines, we use the classical Lick indices in the spectral 
range 4000-5200 A. We introduce a new method for applying resolution corrections 
to the line-strength indices. We define a subset of galaxies with very low emission 
contamination, based on the Ha line, and fit the index-cr relations for this subset. 
The relations show the continuation of the familiar trends for giant galaxies into 
the low luminosity regime, with little change in slope for most indices. Comparing 
the index— a slopes against predictions from single-burst stellar population models, 
we infer the scaling relations of age, total metallicity, [Z/H], and a-element abun- 
dance ratio, [a/Fe]. To reproduce the observed index— a slopes, all three parameters 
must increase significantly with increasing velocity dispersion. Specifically, we recover: 

Ageoc (7 0.52±0.06±0.10 j Z/Ho( (J 0.34±0.04±0.07 j and a/FeQ( (J 0.23±0.04±0.06 (where the 

second error reflects systematic effects) , derived over a decade baseline in velocity dis- 
persion, a — 30 — 300kms _1 . The equivalent slopes for the subset of galaxies with 
a > 100 km s -1 are similar for age and Z/H. For a/Fe, a steeper slope is recovered 
for the high-cr subset, a/Feoc cr °- 36 ± 07 . The recovered age— a relation is shown to be 
consistent with the observed evolution in the giant-to-dwarf galaxy ratio in clusters at 
redshifts z = 0.4 — 0.8. In a companion paper, we determine the age, [Z/H] and [a/Fe] 
for individual galaxies, and investigate in detail the distribution of galaxy properties 
at fixed velocity dispersion. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
clusters: general 



1 INTRODUCTION 

A key challenge in the study of galaxy formation is to explain 
the emergence and build-up over cosmic history of the tight 
red sequence of non-star-forming galaxies, especially in the 
richest environments. At the present epoch, red sequence 
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galaxies account for some 50 per cent of the total stellar 
mass (e.g. Bell et al. 2004). 

The evolution of the red galaxy luminosity function pro- 
vides a direct probe of the mass-assembly history of the red 
sequence. Recent surveys suggest growth in the red sequence 
mass by factors of 2-3 since redshift 0.7 — 1.0, dominated 
by galaxies with relatively low luminosity ( < L*), while the 
most luminous red galaxies show little number evolution 
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over the same interval (e.g. Bell et al. 2004; Faber et al. 
2005; Bundy et al. 2006; Wake et al. 2006; Brown et al. 
2007; Scar lata et al. 2007). The history of star formation is in 
general distinct from the mass-assembly history of the galax- 
ies, since dissipationless ('dry') mergers of gas-poor systems 
merely redistribute existing stars (e.g. van Dokkum 2005). 
Faint field-galaxy surveys reaching to high redshifts reveal 
that star formation has shifted from high-mass galaxies in 
the past to being preferentially active in low-mass systems 
today, a trend referred to as 'downsizing' (e.g. Cowie et al. 
1996; Juneau et al. 2005). 

In galaxy clusters, the Butcher-Oemler effect (Butcher 
& Oemler 1984) constituted some of the earliest evidence for 
the recent rapid growth in the red sequence: a population 
of blue star-forming galaxies existed in some rich clusters at 
z ~ 0.3 which is apparently absent at the present day, having 
presumably faded into the quiescent population in the inter- 
vening ~3 Gyr. Rapid evolution in luminosity-selected stud- 
ies is driven largely by a decline in low-mass star-forming 
galaxies, which fade onto the faint end of the red sequence 
(e.g. Holden et al. 2006). After star-formation is 'quenched', 
especially if the remaining gas was consumed in a rapid star- 
burst, the transformed galaxies may pass through a K+A 
phase (e.g. Dressier & Gunn 1983; Couch & Sharpies 1987). 
The characteristic mass scale of K+A galaxies is observed 
to be larger at high redshift (e.g. Tran et al. 2003) than in 
local clusters (e.g. Poggianti et al. 2004) , a manifestation of 
downsizing in the formation of the red sequence itself. Many 
recent studies suggest that distant clusters are deficient in 
faint red galaxies, relative to local analogues (e.g. Kodama et 
al. 2004; de Lucia et al. 2004, 2007; Muzzin et al. 2007; Stott 
et al. 2007) , as expected if many of today's faint red sequence 
galaxies were indeed star-forming at look-back times above 
a few Gyr (but see Andreon 2006 for contrary evidence). If 
the present day red population is compared naively to the 
galaxies already on the red sequence at high redshift, we in- 
evitably find only modest evolution the stellar populations 
(the so-called 'progenitor bias'), since we are not comparing 
to the full progenitor set (e.g. Kaviraj et al. 2006). None 
the less, recent studies of the red-sequence scaling relations, 
such as the Fundamental Plane, indicate a faster evolution 
in the stellar populations of lower luminosity galaxies, again 
consistent with the downsizing scenario (e.g. J0rgensen et 
al. 2006). 

An alternative means of probing the build-up of the 
quiescent population is to measure explicitly the stellar pop- 
ulation ages of galaxies on the red-sequence at the present 
day. A mass dependence in the quenching times should show 
up in the current ages of galaxies, so long as the latter are 
not in general contaminated by ongoing or renewed star- 
formation significantly after the primary quenching event. 
(This assumption might be better justified in cluster envi- 
ronments than in the low-density field where rejuvenation 
through mergers could be more common.) Determining the 
average stellar ages of galaxies from their integrated light 
is fraught with difficulties, from the classic age-metallicity 
degeneracy (e.g. Worthey 1994) to the presence of non-solar 
abundance ratios (e.g. Worthey, Faber & Gonzalez 1992) 
and the confounding effects of composite stellar popula- 
tions (e.g. Serra & Trager 2007). Despite these challenges, 
modern spectral synthesis models, notably those of Worthey 
(1994) and Thomas and collaborators (Thomas, Maraston 



& Bender 2003; Thomas, Maraston & Korn 2004), have en- 
abled great progress in separating the effects of age, metal- 
licity and abundance ratios on the principal spectral lines. 
(While we will focus here on line-index methods, we note 
that Heavens et al. (2004) have developed an alternative, 
'full-spectrum', method to constrain star-formation histories 
from the present-day galaxy population.) 

Many studies have employed these and similar models 
to analyse modest samples of low-redshift galaxies, mostly 
corresponding to giant elliptical and SO galaxies. We high- 
light here only a selection of notable work. Kuntschner 
(2000) concluded that Fornax cluster ellipticals form a 
metallicity sequence at constant (old) age, while SOs in 
Fornax show a wider age spread. Trager et al. (2000b) 
recovered a wide spread in ages, with an age-metallicity 
anti-correlation maintaining the tight broad-band colour- 
magnitude relation. They show a tendency for younger ages 
at lower a, at least for galaxies outside of rich clusters. More 
recently Thomas et al. (2005) compiled data from a num- 
ber of sources for ~200 galaxies, arguing for a fairly shallow 
age-mass relation (although this conclusion was reached af- 
ter excluding some young low- mass objects from their fits). 

Complementary to these high signal-to-noise (S/N) 
studies of small galaxy samples, a number of works have 
used spectroscopic surveys of thousands of galaxies to ex- 
tract the age-mass relations for red galaxies, although with 
lower typical S/N per galaxy. Among these, Nelan et al. 
(2005), working from the red-sequence-selected NO AO Fun- 
damental Plane Survey (NFPS), strongly confirmed the ear- 
lier hints at an age-mass trend, with Ageoc (j - 6 * 01 ( see 
also Smith et al. 2006). Subsequently, similar trends were 
obtained by selecting red and/or early- type galaxy sam- 
ples from the SDSS: Bernardi et al. (2006) recovered Ageoc 
^.0.8-1.2 ((} e p en(; i m g on w hich age estimator is used), while 
Clemens et al. (2006) found Age oc a 0,8 (our estimate from 
their figure 10). Although these surveys suggest an emerging 
consensus on the existence of the age-mass relation, there 
are still some apparently contradictory results. For instance, 
Sanchez-Blazquez et al. (2006b) report uniformly old ages in 
their sample of ~20 Coma ellipticals (including five galaxies 
with a = 30 — 100 kms -1 ), although their Virgo and 'field' 
samples show a trend with Ageoc a 0,5-10 . It is possible 
that different selection methods, in particular colour selec- 
tion versus detailed visual morphological inspection, are at 
the root of these differences. 

Given the strong intermediate-redshift evidence for 
downsizing and the mass-dependent build-up of the red se- 
quence, there is a pressing need to determine present-day 
galaxy ages over a wide mass baseline. Among the studies 
described above, however, most are limited to galaxies near 
the break in the luminosity function, M* . For example, the 
median velocity dispersion, a, is ~200kms~ 1 for the SDSS- 
based studies (Bernardi et al., Clemens et al.), and for the 
Thomas et al. compilation (which includes only four galax- 
ies with a < 100 kms -1 ). The NFPS is somewhat better in 
this regard, with median a around 150 kms -1 , but the S/N 
achieved on individual galaxies is sufficient only to deter- 
mine an average age at low mass. To trace the build-up of 
the red sequence since z ~ 1, where the activity is mostly 
in faint galaxies, it is necessary to extend the sampling to 
still lower a, while maintaining high spectral S/N (required 
for reliable age measurements), and samples of meaningful 
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size and reproducible selection criteria. A first step towards 
this goal was made by Mobasher et al. (2001) in their sur- 
vey of Coma cluster members, extending to low luminosities. 
However, the resolution of this survey was not sufficient to 
measure velocity dispersion to act as a mass proxy, and the 
S/N (per Angstrom) was too low to determine meaningful 
ages for individual galaxies in the fainter part of the sam- 
ple. Using the Mobasher et al. data, Poggianti et al. (2001) 
found only a very weak dependence of age on luminosity, and 
a very wide spread in age at all luminosities. Caldwell, Rose 
& Concannon (2003) made a higher-resolution and higher- 
S/N study of 175 galaxies (many in Virgo), with 40 galaxies 
in the regime a = 40 — lOOkms -1 , recovering a trend of 
Ageoc cr ' 8 . More recently, Matkovic & Guzman (2005) have 
observed ~90 galaxies in Coma, including 37 galaxies with 
a < lOOkms -1 and S/N > 20 per Angstrom. This sample 
appears to follow an age-mass relation with a similar slope 
to the Nelan et al. result (Matkovic, private communica- 
tion). Finally, two studies of true dwarf galaxies illustrate 
the wide range in properties which perhaps co-exist in the 
low-luminosity population: Geha, Guhathakurta & van der 
Marel (2003) observed 17 Virgo dwarf ellipticals (dEs) with 
o = 23 — 45 kms -1 , and obtained an average age of ~5 Gyr. 
At the other end of the surface-brightness scale, a sample of 
six ultra-compact dwarfs in Virgo, also with a ~ 30 kms" 1 , 
yields ages of > 10 Gyr (Evstigneeva et al. 2007) . 

This paper presents the first results from deep observa- 
tions of galaxies in the Shapley Supercluster with the Anglo- 
Australian Telescope (AAT). The objective is to investigate 
the not only the mean scaling relationships but also the in- 
trinsic scatter in the stellar population parameters, their 
correlation structure, and their dependence on luminosity, 
morphology and environment. Particular emphasis will be 
placed on the faint red-sequence galaxies, which are expected 
to have experienced the most recent truncation events. Such 
an analysis requires age and metallicity measurements for 
hundreds of individual galaxies, spanning a wide range in 
mass, and reaching the level of precision so far achieved only 
for small samples of high-mass objects. The Shapley Super- 
cluster (Proust et al. 2006 and references therein), is an ideal 
target for such a campaign, since it provides an exceptional 
density of galaxies and clusters within a few square degrees 
on the sky. The angular scale of the supercluster and the 
density of feasible targets are well matched to the two de- 
gree diameter field of the AAOmega fibre-fed spectrograph 
(Sharp et al. 2006). For a sample with these characteristics, 
AAOmega delivers a larger multiplex advantage than any 
instrumentation available on larger telescopes. 

The structure of the paper is as follows. In Section[5]we 
describe the sample selection (|2.1[) . the spectroscopic obser- 
vations and data reduction (|2.2[) . and measurements of red- 
shift and velocity dispersion (|2.3p . emission line equivalent 
widths l|2.4p . and absorption line indices (I2.5f) . In Section [3] 
we present statistics of objects showing significant emission 
lines. We next define the sample of low-emission supercluster 
members for further analysis, and tabulate data for this sam- 
ple (Section |4]). Section [5] presents an initial analysis of the 
stellar population scaling relations, based on the index— a 
relations. A comparison of these results to previous work is 
made in Section [6] and the the main conclusions reviewed 
in Section [7] 

In a companion paper (Smith, Lucey & Hudson 2007, 



hereafter Paper II), we use our line-strength data to infer 
the age, total metallicity and [a/Fe] ratio for each galaxy, 
and investigate the age-metallicity-mass relations and their 
intrinsic scatter. 

Throughout this paper, we adopt cosmological param- 
eters (Qm ,^a, h) = (0.3,0.7,0.7). For reference, at the red- 
shift of Shapley (z — 0.048), one arcminute corresponds to 
57kpc, and the distance modulus is m — M — 36.65. 



2 OBSERVATIONS AND PARAMETER 
MEASUREMENTS 

2.1 Sample selection 

The primary galaxy sample was selected using BR imaging 
data from the NOAO Fundamental Plane Survey (NFPS), 
described by Smith et al. (2004). The images were obtained 
with Mosaic-II at the Cerro Tololo Inter- American Obser- 
vatory 4m telescope in March 2000. Targets were selected 
from the same photometric catalogues, generated using SEx- 
tractor (Bertin & Arnouts 1996), as were used to defining 
the NFPS spectroscopic sample. The systematic zero-point 
errors in the catalogues are estimated at ~0.03 mag. Star- 
galaxy separation is via the usual stellarity index of SExtrac- 
tor, limited by image quality ranging from 0.9 to 1.5arcsec 
FWHM. 

The input catalogue for the spectroscopic phase was 
defined purely by total magnitude with R < 18 (one mag- 
nitude deeper than NFPS), before correction for galactic 
extinction (median E(B - V) = 0.053, A R = 0.14, from 
Schlegel, Finkbeiner & Davis 1998). No cuts were made on 
the colour relative to the red sequence (unlike NFPS), nor 
on morphological type. No redshift information was used in 
selecting the sample, so some contamination by background 
and foreground galaxies is expected. The input catalogue 
is inevitably incomplete with respect to low surface bright- 
nesses (SB), since very low-SB galaxies will not have been 
detected by SExtractor. Moreover, very compact galaxies 
will be misidentified as stars in the imaging data available 
for selection (M32, placed at the distance of Shapley, would 
have an effective radius of ~0.1arcsec). Figure \T\ shows the 
target catalogue in a B — R versus Mr colour-magnitude di- 
agram. Pre-empting the following sections, the figure distin- 
guishes the galaxies to which fibres were ultimately assigned, 
and also those which prove to be within the supercluster and 
have pure absorption-line spectra. 

The NFPS imaging data were obtained in 
~40x40arcmin 2 tiles centred on the three rich clusters 
Abell 3556, Abell 3558 and Abell 3562 (for comparison, the 
virial radius of a a c \ = 1000 kms -1 cluster is ~ 45arcmfn 
at the Shapley distance). As such, they cover only a small 
part of AAOmega's two degree diameter field of view. 
Due to fibre-placement limitations, the target galaxies do 
not exhaust the available fibres. Supplementary objects 
were added from the 2MASS Extended Source Catalogue 
(Jarrett et al. 2000), which is limited approximately by 
J < 15, equivalent to _R<15.6 for red galaxies. These 
galaxies are more extended over the field, and could readily 
be assigned unused fibres. However, to maintain a very 
simple sample description, the 2MASS targets will not be 
analysed in this paper. 
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Figure 1. Sample selection from the optical colour magnitude 
relation using NFPS photometry. Colours here are differences of 
total magnitudes, not true matched-aperture colours. The data 
are shown after correction for galactic extinction. The filled circles 
show the galaxies which are used in the absorption-line analysis 
(see Section [3); open circles show galaxies in the redshift range 
adopted for Shapley, but rejected from our sample based on emis- 
sion at He*. Open squares represent foreground and background 
galaxies based on our spectra, while crosses denote objects in 
the NFPS imaging regions but to which spectrograph fibres were 
not assigned. The dashed and dotted lines show the photometric 
catalogue limits of R < 18 and B < 20, for the median galactic 
extinction in the sample, while the arrow shows the characteristic 
luminosity M* for red galaxies in the SDSS (Baldry et al. 2004). 
The solid line shows a robust fit for the emission-free superclus- 
ter members. In the lower panel, we show the fraction of R < 18 
galaxies observed, relative to the NFPS photometric catalogue 
(dashed line), and the fraction of NFPS catalogue galaxies which 
are used in the absorption-line analysis (solid line). 



2.2 AAOmega spectroscopy 

Spectroscopic observations were made with the AAOmega 
system (Sharp et al. 2006) on the 3.9m Anglo- Australian 
Telescope, on the nights of 2006 April 26 and 29. AAOmega 
consists of the 2dF fibre positioner (Lewis et al. 2002), linked 
to an efficient and stable bench-mounted dual-beam spectro- 
graph. The positioner can place 392 fibres, with 2 arcsec pro- 
jected diameter, within a ~2 degree diameter field of view. 
Within each configuration there is a minimum target sep- 
aration of ~30 arcsec imposed by the physical size of the 
fibre buttons (for further details of the fibre positioning con- 
straints, see Miszalski et al. 2006). 



Two fibre configurations were employed, with substan- 
tially overlapping fields of view. The allocation of fibres to 
the targets was performed in a two-stage process. First, ap- 
proximately 70 fibres known to exhibit significant fringing 
effects were withdrawn, and the remaining 'good' fibres al- 
located to the 72 < 18 sample galaxies. Some 30 'good' fibres 
were assigned to random positions to sample the background 
sky spectrum. In a second step, the locations of the these 
fibres were held fixed, while remaining fibres (both 'good' 
and 'fringed') were allocated where possible to extra targets 
from the 2MASS sample. 

In the blue arm of the spectrograph, the 580V grating 
was used, delivering a nominal resolution of 3.5A FWHM 
from 3700-5800 A, sampled at TO A per pixel. The blue data 
thus cover most of the standard Lick absorption lines, and 
also the Oil 3727A emission line when present. Simultaneous 
observations were obtained in the red arm of the spectro- 
graph, using the 1000R grating, in order to measure nebular 
emission at Ha. The instantaneous wavelength coverage is 
~1100A and the nominal resolution is 1.9 A, sampled at 
0.6 A per pixel. The two fields were observed alternately, in 
visits of 1.5-3.0 hours, to a total integration time of 8-9 
hours. In this way, each configuration was observed always 
with the same field plate and fibre bundle, so that all the 
spectra for a given galaxy were obtained with the same fi- 
bre. (Exceptions are seven objects common to both fields, of 
which one is a star). Between the two nights, we introduced 
a small shift (~30 A) in central wavelength for the blue spec- 
tra, in order to fill gaps left by blocks of bad columns in the 
detector. In the red arm, a larger shift (~350 A) was used 
to expand the total wavelength coverage (to 5800-7300 A). 
A total of 565 targets were observed, with 416 from the pri- 
mary NFPS sample and 149 from the supplementary 2MASS 
list. Overall, ~60% of the NFPS sample targets were ob- 
served. Note that because the magnitude information was 
not used in prioritising fibre assignment, the observed frac- 
tion is not a strong function of luminosity. 

The data were reduced using a combination of the 
AAO drcontrol pipeline with custom-written procedures. 
The two-dimensional frames were first cleaned of manually- 
identified bad columns, while remaining bad pixels and cos- 
mic ray hits were identified using LA COSMIC (van Dokkum 
2001). The well-sampled PSF of AAOmega ensures that cos- 
mic rays are clearly distinguished from narrow sky lines and 
emission lines. The spectra were flat-fielded and extracted 
(using a Gaussian profile extraction) with drcontrol. The 
combination of spectra from the many 30 minute exposures 
is automated, but handled on a galaxy-by-galaxy basis. Ex- 
amination of the data revealed sizable variations in signal- 
to- noise ratio (S/N) between exposures, for some (but not 
all) objects. To lessen the impact of these variations, we first 
determine the maximum S/N attained for the galaxy, then 
discard any individual exposures with S/N less than a quar- 
ter of this maximum value before combining to form the final 
spectrum. This criterion rejects only ~2.5% of the ~10000 
individual spectra obtained, but ~10% of galaxies have at 
least one rejected spectrum. Around half of the affected spec- 
tra are from fibres close to the edge of the field-of-view, 
while the other half are concentrated in an annulus around 
0.27 degrees from the field centre. The S/N variations near 
the edge of the field are more frequent at high airmass, and 
probably result from varying atmospheric refraction over the 
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course of a field visit, which can cause relative positional 
offsets of ~larcsec in the worst cases. For comparison the 
sample galaxies have effective radii typically 2-4arcsec (F. 
La Barbera, private communication). When combining the 
spectra we propagate the error spectra from drcontrol, 
including missing data flags, to produce a final error spec- 
trum for each galaxy. Finally, an approximate relative flux 
calibration was imposed using the observed spectrum of a 
calibration star. 

Each of the combined spectra was examined by eye to 
identify any objects unsuitable for the subsequent analysis. 
These are: seven stars, three objects in which the target is 
clearly contaminated by a second galaxy with very differ- 
ent redshift, or by a star, and 15 galaxies whose combined 
spectra are affected by sinusoidal continuum variations, at- 
tributed to 'fringing'. These objects are rejected entirely 
from further discussion. For the remaining 540 'uncorrupted' 
galaxies, the visual inspection confirms there is sufficient 
signal to measure at least a reliable redshift, and that any 
remaining blemishes (usually from very low-level detector 
defects and visible only in the faintest galaxies) affect only 
localised wavelength regions. Among these galaxies, 402 are 
from the primary (NFPS) sample and 138 from the 2MASS 
supplementary sample. 

Total integrations range from five hours to nine hours 
for most galaxies. Around thirty galaxies have integrations 
from two to five hours, after the S/N cut described above. 
Finally, six galaxies, observed in both field configurations, 
have total integration times of 14-17 hours. 

2.3 Redshifts and velocity dispersions 

Reliable velocity dispersion measurements depend on com- 
paring the galaxy spectra to templates which match them 
as closely as possible. Traditionally, a small number of K 
giant stars have been observed through the same instrumen- 
tal set-up for use as templates. While the red giant branch 
indeed contributes substantially to the integrated V-band 
light of red-sequence galaxies, there is a comparable contri- 
bution from main-sequence stars near the turn-off (e.g. Fig. 
13 of Maraston 2005). The contribution from hotter stars 
varies systematically with wavelength, and also as a func- 
tion of the age and metallicity of the target galaxy itself. As 
a result, a single template star cannot match well over an ex- 
tended wavelength region. Errors associated with this 'tem- 
plate mismatch' can cause systematic bias in the recovered 
velocity dispersions. Methods have been developed which 
minimize template mismatch by fitting optimal combina- 
tions of template stars (e.g. pPXF of Cappellari & Emsellem 
2004). Here however, we adopt an alternative approach, us- 
ing as templates a set of synthetic SSP spectra covering a 
range in metallicity and age. The key requirement is that 
the spectral resolution of the models exceed the instrumen- 
tal resolution. In this paper we employ the synthesis models 
by Vazdekis et al. (2007), based on the 2.4 A MILES spectral 
library (Sanchez-Blazquez et al. 2006c). 

Redshifts are required in order to select a common 
rest-frame wavelength interval for the dispersion measure- 
ments. We determine absorption-line redshifts using the 
cross-correlation method (Tonry & Davis 1981), as imple- 
mented in iraf as FXCOR. The whole observed spectrum is 
used at this stage. A high-order spline fit to the spectrum is 



used to identify strong emission lines, which are interpolated 
over to reduce their influence in the cross-correlation. For ra- 
dial velocity templates, we employ a reduced set of MILES 
SSP spectral models. Specifically, we used solar-metallicity 
SSP models with ages 0.5, 1.0, 2.0, 4.0, 8.0 and 16.0 Gyr, 
adopting the redshift corresponding to the best-matching 
SSP, as judged from the fxcor R-value. The young SSP 
models proved useful in obtaining good absorption-line red- 
shifts from galaxies with substantial emission plus strong 
underlying Balmer absorption. We manually inspected all 
spectra for which the best R-value was less than 10, and 
those with substantial variation in the redshifts derived 
from young and old templates. In all cases, the measured 
redshifts were confirmed. Comparing to redshifts from the 
NFPS (Smith et al. 2004), for 124 galaxies in common, we 
obtain a median offset of 11 kms -1 (the new measurements 
larger), with a scatter of 23 kms -1 . 

As a preliminary to measuring velocity dispersions, we 
determine the delivered resolution of our spectra as a func- 
tion of wavelength, by comparing the spectra obtained for 
a K2 giant, HR7149, with the spectrum for this star in the 
Indo-US stellar library (Valdes et al. 2004). For wavelengths 
of 4500-5750 A, we find a constant 3.0 A FWHM smoothing 
must be applied to the Indo-US spectrum for optimal match 
to the AAOmega spectrum. Adding in quadrature the Indo- 
US resolution (described as 1.2 A FWHM with 'small vari- 
ations'), the AAOmega resolution is thus estimated to be 
3.2 A FWHM (equivalent to <T inst = 82 kms -1 ). For 4000- 
4500 A, the resolution is slightly poorer, ~3.7A. Consider- 
ing the range where the resolution is approximately con- 
stant, and the redshift range spanned by our principal tar- 
gets in Shapley, we select the spectral range 4400-5400 A 
in the rest frame, for determination of the velocity dis- 
persion. This range is fully covered for 461 galaxies with 
cz < 22000 kms -1 . For velocity dispersion templates, we 
again employ the MILES SSP models. In the wavelength 
range of interest, the MILES resolution is 2.4 A FWHM (fig- 
ure 4 of Sanchez-Blazquez et al. 2006c). We must therefore 
smooth the MILES template spectra by the quadrature dif- 
ference (2.1 A), to approximate a zero-velocity SSP as ob- 
served by AAOmega. 

Dispersions were measured using FXCOR again, calibrat- 
ing the cross-correlation peak widths to a using results from 
artificially broadened templates (e.g. see Wegner et al. 1999). 
The trial SSPs form a grid with metallicities [Fe/H]= - 
0.68, -0.38, 0.00, +0.20, and ages 1.0, 1.4, 2.0, 2.8, 4.0, 
5.6, 8.0, 11.2, 16.0 and 17.8 Gyr. As for the redshifts, we 
select the result corresponding to the best-matching SSP 
template. Examination of the results for different templates 
shows that if we had used a single old, solar-metallicity SSP 
template, the dispersions would be strongly biased relative 
to the 'best' SSP, especially for low-mass galaxies. Specifi- 
cally, while for a > 100 kms -1 the two estimates agree with 
an offset of only ~2 per cent, and a scatter of ^11 per cent, 
for a = 75 kms -1 , the median offset is 6 per cent, and for 
a — 50 kms -1 , the bias is 18 per cent. In these cases the 
'best' a is on average smaller than the a obtained from a 
single old SSP template. 

Uncertainties in the velocity dispersions were estimated 
using Monte Carlo simulations using the error spectra cor- 
responding to each galaxy. From 51 realisations, we numer- 
ically sort the resulting measurements and adopt the 8th 
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Figure 2. Signal-to- noise as a function of velocity dispersion for 
the dataset presented in this paper (filled circles). We show only 
the ~300 galaxies with very low Ha emission (see later) . For com- 
parison, we show the sample of Matkovic & Guzman (2005) as 
small crosses, and the NFPS (Smith et al. 2004) as grey points. 
The kinematically unresolved galaxies (i.e. those with velocity 
dispersions consistent with zero) are placed at an arbitrary loca- 
tion at the left of the panel, separated from the other data by a 
vertical line. 

and 42nd values as an estimate of the la interval. In this 
process the best-matching SSP model is used throughout; er- 
rors associated with selecting the template are therefore not 
propagated into the velocity dispersion error. If the galaxy- 
versus-model cross-correlation peak width is narrower than 
the model autocorrelation peak width, at the la lower limit, 
we classify the galaxy as kinematically unresolved. This is 
true for 90 of the total 461 galaxies for which measurements 
were made. The unresolved galaxies are primarily those with 
best match to very young (< 2 Gyr) SSP models; many also 
have strong nebular emission lines. This probably results 
largely from the increased influence of broad A-star Balmer 
absorption, which degrades our effective instrumental reso- 
lution for ages younger than ~3 Gyr. 

Figure [2] shows the relationship between measured ve- 
locity dispersion and the signal-to-noise ratio (S/N, mea- 
sured over the 4400-5400 A rest-frame interval) . The figure 
suggests that the lowest measurable velocity dispersions in 
our data are 30 — SOkms^ 1 (roughly half the instrumental 
resolution of ~90kms _1 ), depending on the S/N achieved. 
For comparison, we show the same relation for the entire 
NFPS spectroscopic dataset (from Smith et al. 2004). For 
galaxies in common between the samples, the new data have 
a factor of ~3.2 higher S/N (~10 times deeper). A more 
comparable sample is that compiled for Coma cluster galax- 
ies by Matkovic & Guzman (2005). Figure [2] shows that at 
comparable a, the S/N in our spectra is larger by a factor 
of around two. 

In Figure [3] we compare the new velocity dispersion 
measurements with those from NFPS (Smith et al. 2004), 



1.4 1.6 1.8 2.0 2.2 2.4 2.6 

log a (NFPS) 

Figure 3. Comparison of velocity dispersion measurements with 
those from NFPS (Smith et al. 2004). The weighted fit to all 
points is shown by the dotted line. Filled points have best- 
matching templates older than 4 Gyr; the dashed line is a fit to 
these galaxies only. The solid line indicates equality. 

for 105 galaxies in common. No aperture corrections need to 
be applied, since the NFPS fibres are of the same diameter 
as the AAOmega fibres. The comparison is limited by NFPS 
to fairly high luminosity. Moreover there is an explicit bias 
because low-a galaxies are likely to be unresolved in NFPS, 
and only enter the comparison if errors cause them to scatter 
up to higher a. We find a median offset of 0.045 dex (~11 per 
cent), with the new measurements consistently smaller. The 
scatter around this offset is 0.068 dex. For low-dispersion 
galaxies, both the scatter and the average offset increase. 
The slope of the (error-weighted) fit is 1.06±0.03, indicating 
a significant scale change. We attribute this behaviour partly 
to the improved measurement method which matches the 
template to each galaxy, as compared to using a few K- 
giant templates as in NFPS. In particular, many of the low- 
er galaxies are best represented by young SSP models. As 
noted above, forcing a fit with a K-giant or old SSP will bias 
the results to higher dispersions in such cases. As a test, we 
repeat the comparison for only those galaxies with a best- 
matching template age greater than 4 Gyr. The agreement 
is considerably improved, with a scatter of 0.048 dex and 
median offset 0.037 dex (~ 9 per cent). Moreover the scales 
are consistent (slope 1.00±0.03). 

Hereafter in this paper, we prefer to use the velocity 
dispersions derived from the best-matching template. A con- 
sequence is that the a scale is slightly stretched relative to 
the NFPS scale, so that measured correlations with a are 
expected to be flatter than if we had only used K-giant tem- 
plates. 

2.4 Emission line equivalent widths 

Since our sample has not been selected by morphology or 
colour, it will contain star-forming galaxies as well as pas- 
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sive systems. In addition, it has long been known that a 
substantial fraction even of morphologically-selected ellipti- 
cals exhibit nebular emission lines such as H/3, [Om] 5007, 
etc (e.g. Phillips et al. 1986). The line ratios may be char- 
acteristic of LINER excitation, especially in more massive 
galaxies, or of excitation by ongoing star formation, the lat- 
ter more frequently in low-mass galaxies (Nelan et al. 2005; 
Yan et al. 2005). 

Emission line measurements are necessary to identify 
passive, non-star-forming galaxies, with minimal nebular 
contamination in the age-sensitive Balmer lines. As dis- 
cussed in Nelan et al. (2005) and elsewhere, the intrinsic 
variation among galaxies in the line ratios precludes accu- 
rate correction for emission infilling based on the [Om] 5007 
line alone, especially if a wide range in galaxy mass is being 
considered. More reliable corrections can, however, be de- 
rived if a wider range of lines are measured, especially if Ha 
is available (e.g. Caldwell et al. 2003). 

To measure equivalent widths of the principal emission 
lines, it is necessary first to remove the underlying stellar 
continuum, especially in the case of the Balmer lines them- 
selves. To this end, we Doppler-corrected each spectrum to 
the rest frame and divided by the best-fitting MILES SSP 
(as determined from the cross-correlation results), broad- 
ened to the measured velocity dispersion. Low-order con- 
tinuum differences were removed using a cubic spline fit. 
Equivalent widths were measured on the ratio spectrum, 
by direct flux summation over an interval equivalent to 
600 km s -1 width centred at the wavelength expected from 
the absorption- line redshift. This measurement window is 
wide enough that modest velocity shifts (— lOOkms ) be- 
tween gas and stars introduce only small variations (—10 
per cent) in the equivalent widths. In Section [3] we will 
use mainly the ratios of neighbouring lines, where the ef- 
fect of such offsets will cancel out. Errors in the equivalent 
widths were estimated from the pixel-to-pixel variance in 
the (template-divided) continuum regions. 

In the blue, the emission lines considered in this paper 
are [On] 3727, H/3, and [Om] 5007, while in the red spec- 
tra, we use Ha and [Nil] 6583. Note that the separation of 
the latter two lines is equivalent to 900kms _1 , i.e. larger 
than the velocity interval used for measuring the equiva- 
lent widths. At the redshift of Shapley, Ha and [Nil] are af- 
fected by atmospheric absorption bands, especially that due 
to 2 at 6870-6955 A. Thus we need to disentangle three 
contributions: stellar continuum, nebular emission and at- 
mospheric absorption. To determine the latter, we use the 
objects with the highest S/N, and with large velocity dis- 
persion, for which any high-frequency variations must be tel- 
luric in origin. The best-fitting SSP models for these galaxies 
were Doppler corrected to the observed frame (the opposite 
of the sense used for continuum removal), broadened to the 
appropriate velocity dispersion and divided out of the spec- 
tra. The ratio spectra for each galaxy retain the atmospheric 
features (always at the same wavelengths) and any emis- 
sion lines present (at different wavelengths according to the 
redshift) . Median-combining the ratio spectra effectively re- 
moves the emission component, leaving a smooth continuum 
modulated by the telluric absorption. The final step was to 
fit out the continuum and create a correction spectrum fixed 
to unity for all pixels below 6800 A. (We have neglected the 
much weaker absorption features at shorter wavelengths). 



Telluric corrections were derived independently for the two 
fields, since these cover different time intervals. A compari- 
son shows that both the overall amplitude and the detailed 
structure of the O2 band (and the H2O band longwards of 
7130 A) are very similar between the two fields. 

2.5 Absorption indices 

The Lick absorption-line indices from HdA0 to Fe5406 were 
measured using indexf (see Cenarro et al. 2001), on the 
flux-calibrated spectra, at the native spectral resolution of 
-3.2 A. The indices Mgl, Mg2, Fe5270 and Fe5335 are not 
reported here, since at the redshift of Shapley, these are con- 
taminated by the 5577 A sky emission for a large fraction of 
the galaxies. We use the passband definitions of Trager et 
al. (1998), and the the conventional units, i.e. magnitudes 
for CN1, CN2 and angstroms for all others. Uncertainties 
are estimated from the error spectrum for each galaxy. Be- 
fore the indices can be compared against spectral synthesis 
models, three distinct corrections must be considered: cor- 
rection to zero-velocity broadening, correction to the spec- 
tral resolution of the models, and correction to match the 
instrumental response curve used in the models, if this is 
not a fluxed system (e.g. the Lick/IDS system). 

2.5.1 Velocity broadening corrections 

The measured absorption line indices differ from those com- 
puted by spectral synthesis models, due to broadening by the 
line-of-sight velocity distribution function. For an isolated 
absorption line, increasing a causes the line to 'leak' from 
the central index band, depressing the measured index. With 
further broadening, the absorption itself spreads into the 
pseudo-continuum bands, further reducing the measured ab- 
sorption. In a complex spectrum with many overlapping fea- 
tures, the effect of velocity broadening can be assessed only 
empirically by comparison to artificially-broadened template 
spectra. 

Most studies in the field have derived velocity broad- 
ening corrections (VBCs) to the indices, based on measure- 
ments on artificially broadened spectra of template stars ob- 
tained with the same instrumentation as the galaxies (e.g. 
Gonzalez 1993; Davies, Sadler & Peletier 1993; Trager et al. 
1998). Alternatively, a library of SSP template models can 
be used to cover a wider range in spectral type, and provide 
a better match to the composite spectra of real galaxies 
(e.g. Kuntschner 2004). The correction described by each 
template spectrum is usually normalised by the index value 
at zero velocity dispersion, and the results from different 
templates averaged into a single multiplicative correction 
describing the ratio A number of indices (most notably the 
H7 and H<5 indices of Worthey & Ottoviani 1997) can take 
a range of values which span or approach zero. Under such 
conditions the usual multiplicative correction fails, and in 
these cases an additive correction, based on 1(a) — I(o- — 0) 
has generally been used. 

In general, it is not clear whether one should apply an 

1 We use roman type, e.g. HdA for the names of the conventional 
line-strength indices, while H<5 refers to the line itself, whether in 
absorption or emission. 
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additive or a multiplicative correction for any given index. In 
the simplest case, where a single correction curve is used for 
all the target galaxies, a multiplicative correction may intro- 
duce systematic errors as a function of linestrength (and thus 
indirectly as a function of a). More recently, many works 
have derived corrections for each galaxy based on a model 
template chosen to match that galaxy (e.g. Kuntschner et al. 
2006; Sanchez-Blazquez et al. 2006a; Kelson et al. 2006). In 
this paper, we adopt an alternative approach to the VBCs, 
which allows a general and consistent treatment of all the 
indices. 

Our corrections are based on a range of MILES SSP 
model templates (Vazdekis et al. 2007), with ages 3-18 Gyr 
and metallicities [Fe/H] = -0.38, 0.00 and +0.20. We smooth 
the template spectra to mimic both the AAOmega instru- 
mental resolution, and also the degradation caused by ve- 
locity broadening in steps from zero to 500 km s" 1 . (We ne- 
glect any higher order moments of the line-of-sight velocity 
distribution, but in principle the method could be general- 
ized to allow for these.) The indices are measured on the 
broadened spectra. For each step in a, the broad coverage 
of spectral types yields a wide range in index value. At each 
value of u, we then simply compare the index measured on 
the velocity broadened spectrum, lNat(c), with that mea- 
sured only allowing instrumental resolution, iNat(c = 0). 
(The subscript 'Nat' here refers to measurements at the na- 
tive resolution of the spectra.) For most of the standard 
indices, this comparison yields a closely linear relation such 
that /Nat(cr = 0) = o/n at (cr) + b (at fixed a), with very small 
scatter. This can be applied as a correction for the indices 
measured on the galaxy spectra. The coefficients a and b for 
each galaxy can be interpolated from the steps in a used for 
smoothing the templates. In most cases the SSP models span 
the range of the measured values. Where this is not the case 
(e.g. the highest Mgb5177 values, for a-enhanced galaxies), 
the high degree of linearity in the corrections suggests that 
a small extrapolation is probably safe. In applying the cor- 
rections, we account for the effect of the linear correction to 
the index error. The uncertainty in the velocity dispersion 
(and hence the coefficients a and b) could also be propa- 
gated to the index error, but in practice this contribution is 
negligible compared to the index error. 

The left panel of Figure [3] shows the corrections derived 
for the Mgb5177 index, for a = 100, 200, 300 kms~\ In this 
case, the correction lNat(o" = 0) — /Nat(f") increases with 
index strength such that our correction is almost identical 
to the usual multiplicative form. In Figure [5] we present the 
best linear corrections at a = 300 km s" 1 , for a selection 
of the principal indices, with a comparison to the equiva- 
lent pure multiplicative and pure additive corrections. The 
figures demonstrate that the new approach yields linear cor- 
rections which hold to a high degree of precision over a very 
wide range of spectral types. By contrast, for some indices, 
the best multiplicative and best additive corrections can be 
accurate only over a limited range of the index value itself. 

As an example, consider the fairly narrow iron index 
Fe5270. The linear and multiplicative corrections (as usu- 
ally employed for this index) reproduce the observed values 
over the whole range of SSPs. By contrast, in the panel for 
HdF, we find that the multiplicative correction fails badly, 
as expected for this index. The additive correction, as usu- 
ally employed, is as successful as our linear fit, although 



in both cases there is some metallicity dependence in the 
residuals. In these cases, our results support the form of 
correction generally adopted. Two counter-examples, where 
our scheme suggests a different form of correction than the 
standard method, are HgF and G4300. For HgF, the tra- 
ditional additive correction is a poor fit to the broadened 
spectra. Our linear correction is much closer to the mul- 
tiplicative form, but is a better fit than multiplicative for 
the most Balmer-strong (young or metal-poor) models. For 
G4300, previous work has generally applied a multiplicative 
correction. Our analysis shows that an additive correction 
is much preferred for this index. However, the index has 
a complicated metallicity-dependent behaviour which is not 
accurately reproduced even by our linear correction scheme. 

In summary, the new formulation of the VBCs described 
here, while not much more complicated or difficult to im- 
plement than the traditional approach, improves upon the 
traditional method in at least two regards: (i) greater gen- 
erality, requiring no subjective preference for additive or 
multiplicative corrections, and (ii) ready applicability over 
a wide range of spectral types for most indices. Our ap- 
proach should be equivalent to the more modern methods 
which derive corrections from matched template spectra. We 
contend, however, that describing the VBCs in terms of a 
mapping from observed index to the value at a — is a more 
natural way to present the correction. 

2.5.2 Instrumental resolution correction 

In order to compare to models defined at the resolution of 
the Lick/IDS stellar library, such as those of Worthey (1994), 
Thomas et al. (2003, 2004), and Schiavon (2007), it is neces- 
sary to match the 8-12 A wavelength-dependent spectral res- 
olution of the IDS instrument. Traditionally this is achieved 
by artificially broadening the spectra prior to measuring the 
indices. A serious drawback is that this results in redistribu- 
tion of any non-uniform noise in the spectra. For example, in 
intermediate-redshift observations, residuals from subtract- 
ing narrow sky lines would be smeared over ~20 A intervals. 
In the data described here, noise non-uniformity is less se- 
vere, but arises from the blocks of bad columns in the detec- 
tor. Since our observed wavelength range was shifted mid- 
way through the observations there are (usually) no pixels 
with unknown values, but in certain regions the variance is 
doubled since only half the integration yields good data. 

An alternative approach, as pointed out by Kelson et 
al. (2006), is to subsume the resolution correction into the 
same process as the VBCs. The correction to the Lick res- 
olution can be obtained by measuring indices on the tem- 
plate SSP models after broadening to the Lick resolution, 
but without any velocity broadening, denoted /Lick(c = 0). 
The smoothing required for each index is as tabulated by 
Schiavon (2007), based on fig. 7 of Worthey & Ottoviani 
(1997). Corrections from observed to zero velocity disper- 
sion can then be established at the Lick instrumental reso- 
lution, in the same way as before, by fitting a linear relation 
^Lick(o" = 0) = a/Nat (c) + b, at fixed a The right panel of 
Figure [4] shows the correction for Mgb5177 as an example. 
Since the Lick resolution is equivalent to a ~ 200kms _1 , 
our higher-resolution measurements are subject to a nega- 
tive correction at low a and to a positive correction at high 

<T. 



© 2007 RAS, MNRAS 000. [114201 



Shapley Supercluster: Stellar population trends 9 



Mgb5177 correction to sigma=0 at native resolution 



Mgb51 77 correction to sigma=0 at lick resolution 



| i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i 




100 

I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Observed lNat(°) 



| i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | 



Fe/H 
Fe/H 
Fe/H 



■■ -0.38 

: 0.00 
: +0.20 



300 




200 



100 

I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ' ■ I ' 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Observed lNat(a) 



Figure 4. Left panel: velocity broadening correction for Mgb5177, correcting to a = at the native resolution of the spectra. For 
each index, we plot /Nat( CT = 0) — /Nat( CT )i where iKat( CT = 0) i s the index measured on the MILES model SSP spectra with no 
velocity broadening (but smoothed to our native instrumental resolution), and /Nat(c) the index measured after broadening to a = 
100, 200, 300 km s _1 , as well as to the native AAOmega resolution. The latter is thus the index value as 'observed'. Points show the 
corrections for SSPs with [Fe/H]=— 0.38 (squares), 0.00 (triangles) and +0.20 (circles). The thick grey line shows zero correction. Solid 
lines are linear fits to the relation given by the MILES SSPs. Right panel: Equivalent corrections to a = at the Lick resolution. In this 
case, /Nat( cr ) is, as before, the index after broadening by a the native instrumental kernel, while /Lick( cr = 0) is measured with the Lick 
instrumental smoothing, but with no velocity broadening. 



To recap, our corrections to zero velocity dispersion and 
to the Lick resolution are handled in a single step. The 
corrections can be applied without smoothing the observed 
galaxy spectra, correctly preserving the noise properties of 
the data. 



NFPS, which are typically three times larger than those in 
the new measurements. In the case of Mgb5177, which is 
well measured in both studies, the data show a zero-point 
offset of 0.10±0.04. In the other indices shown, no significant 
offset is observed. 



2.5.3 No calibration to the Lick system 

Finally, the most widely used models of the past few years 
are those of Thomas et al. (2003, 2004), which are based 
on the Lick/IDS stellar library (Worthey et al. 1994). This 
library is not flux calibrated, instead being defined by the 
wavelength response of an instrument of only historical sig- 
nificance. We did not expend valuable telescope time observ- 
ing Lick standard stars to match this system. Instead, when 
we compare to the Thomas et al. models we will employ 
them in a relative sense only. This can be done either by 
considering only the slope of, for instance, the index— a re- 
lations (as in Nelan et al. 2005), or by imposing zero-point 
shifts to fix a set of the most massive galaxies to match 
fiducial values for age, [Z/H] and [a/Fe]. Ultimately, we an- 
ticipate re-analysing the current data in the context of a 
new generation of models now being developed (e.g. Schi- 
avon 2007; Vazdekis et al. 2007) which are tied to carefully 
flux-calibrated stellar libraries (e.g. Jones 1999; Valdes et al. 
2004; Sanchez-Blazquez et al. 2006c). 

2.5.4 Comparison with NFPS 

Figure HJ] compares the new index data with that of NFPS 
(Nelan et al. 2005) for 78 galaxies common to the two sam- 
ples. The comparisons are limited by the uncertainties in 



3 EMISSION LINE ANALYSIS 

In this section, we analyse the emission line equivalent 
widths, primarily to determine a criterion for selecting an 
emission-free galaxy sample. The incidence of emission lines 
in cluster members is also of interest in itself, reflecting the 
ongoing consumption of gas in the cluster environment. 

In the following sections, we will use indices measuring 
the Balmer lines (specifically Hbeta, HgF, HdF) as age indi- 
cators for the galaxy sample, on the assumption that these 
are dominated by absorption in stars near the main-sequence 
turn-off. If the galaxy contains regions of star formation, or 
hosts an active galactic nucleus (AGN or LINER), then the 
Balmer lines are contaminated by emission, yielding weaker 
net absorption and spuriously old age estimates. Such galax- 
ies must either be corrected for the effects of emission or 
excluded from the sample entirely. Many older studies chose 
to correct for H/3 emission by assuming a constant ratio 
between the [Oni] 5007 and H/3 lines (e.g. H/3/[Om]= 0.6, 
Trager et al. 2000a). However, the measured line ratios, 
even for red-sequence galaxies, span a wide range around 
this value, and moreover vary systematically with galaxy 
mass (Nelan et al. 2005; Yan et al. 2006). Specifically, cor- 
rection using [Om] under-corrects the emission at low mass 
and over-corrects at high mass. Using [Om] alone to select 



© 2007 RAS, MNRAS 000. 01-1201 



10 Russell J. Smith et al. 



HdF rms = 0.018 rms(x) = 0.169 rms(+) = 0.019 



i ■ ■ ■ ■ i ■ ■ 11 i ■ ■ ■ ■ i 11 11 i ■ ■ 1 1 1 1 ; 



MA A m □ 
-g O aa^fflAAAoA^;^" j o UP""— 



I .... I I , , , , I , , , , I , , , , I 



-1.0 -0.5 0.0 0.5 1.0 1.5 
Observed lNat(") 



Ca4227 rms = 0.012 rms(x) = 0.022 rms(+) = 0.139 
1 



G4300 rms = 0.026 rms(x) = 0.077 rms(+) = 0.034 

r 




1.0 1.5 2.0 

Observed l Na t(<i) 



4.0 4.5 5.0 5.5 
Observed iNatM 



Fe4383 rms = 0.010 rms(x) = 0.068 rms(+) = 0.155 

r 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 
Observed lNat(c) 



4 5 
Observed l Na t(o) 



Hbeta rms = 0.020 rms(x) = 0.020 rms(+) = 0.029 
T 




1.5 2.0 2.5 

Observed iNat(a) 



Fe501 5 rms = 0.026 rms(x) = 0.040 rms(+) = 0.207 



Mgb5177 rms = 0.008 rms(x) = 0.011 rms(+) = 0.094 



Fe5270 rms = 0.007 rms(x) = 0.015 rms(+) = 0.110 

r 




4 5 
Observed lNat(a) 



2.5 3.0 3.5 4.0 
Observed 



2.0 2.5 3.0 3.5 
Observed lNat( n ) 



Figure 5. Velocity broadening corrections, to the native AAOmega resolution, for a selection of important indices, as derived for 
a = 300kms — 1 . The thick grey line shows zero correction. The solid black line is a linear fit to the relation given by the MILES SSPs, 
i.e. the correction used in this paper. The dashed line (often very nearly coincident with the solid) shows the best-fitting multiplicative 
correction, i.e. a constant factor I (a = 0)/I(a = 300) for all spectral types. The dotted line shows the best fitting-additive correction, 
i.e. constant offset I(cr = 0) — I (a = 300) for all spectral types. In the title bar, we note the RMS of the models around the best linear 
relation, as well as the RMS around the multiplicative and additive corrections. 



low-emission samples risks the inclusion of high H/3/[Oni] 
objects at low a. Either case could result in a spuriously 
flattened age— a relation. 

A more useful diagnostic, when available, is a measure- 
ment of Ha emission, since the Ha/H/3 ratio should be much 
more stable than, for instance [Oni]/H/3, which depends 
strongly on the source of excitation. Moreover, nebular emis- 



sion at Ha dominates over the stellar absorption even for 
very low-levels of star formation or AGN luminosity, provid- 
ing powerful leverage on the contamination of higher-order 
lines, and robustness against errors in continuum removal 
(e.g. Caldwell et al. 2003). 

The sample used in this section comprises 342 galaxies 
selected from the NFPS imaging, with redshifts compatible 
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Figure 6. Comparison of the new index measurements with data from NFPS (Nelan et al. 2005), for three representative indices. There 
are 78 galaxies in common which meet our emission criterion (EW(Ho)<0.5 A). The line indicates equality. 



with membership of the supercluster (velocities in the range 
11670-17233 km s _1 , justified in the following section). Ha 
emission is detected with EW(Ha)> 0.5 A, for 110 objects 
i.e. 32 per cent. With this threshold, the detection of Ha is 
always significant at the > 3<r level. Of these 110, emission 
is also detected with S/N>3 in H/3 for 85 per cent, in [Oni] 
5007 for 89 per cent, in [Nil] 6583 for 96 per cent and in [On] 
3727 for ~80 per cent. The converse proportions are more 
critical: we are most concerned to ensure that the measured 
H/3 absorption is not substantially contaminated by nebu- 
lar emission. We prefer not to use the directly-measured H/3 
emission to impose a sample cut, so that errors in selecting 
the best-fitting continuum template are not propagated di- 
rectly to the absorption line measurements. However, we can 
use these data to assess the reliability of selecting accord- 
ing to the other lines. For 232 galaxies without detected 
Ha ('Ha-clean'), only ten (i.e. 4 per cent ) show signifi- 
cant (S/N>3) H/3 emission. The equivalent fraction for 269 
[Om]-clean objects is 10 per cent, while for 241 [Oll]-clean 
objects, 12 per cent show H/3 emission. (In defining clean 
samples from [On] and [Oni], we did not impose an abso- 
lute EW cut comparable to the 0.5 A applied for Ha, since in 
these cases there is not a strongly- varying underlying stellar 
absorption line.) Given the direct astrophysical connection 
between Ha and H/3, we would have expected this to be 
the most accurate indicator for nebular contamination. The 
tests presented in this paragraph demonstrate in addition 
that it is also the most sensitive indicator from a practical 
perspective, admitting a factor 2-3 times fewer galaxies with 
H/3 emission than the available lines in the blue spectra. 

Next we consider the origin of the emission lines, i.e 
whether they are due to star formation or to AGN/LINERs. 
The upper left panel of Figure[7]shows the classical emission- 
line diagnostic plane of Baldwin, Phillips & Terlevich (1981), 
hereafter the BPT diagram. (Note that while ratios of line 
fluxes are usually used in such diagrams, we use ratios of 
equivalent widths. Since we restrict this analysis to line- 
pairs which are close in wavelength, the ratios differ only at 
the ~ 0.02 dex level from flux ratios.) Normal star-forming 
galaxies are distinguished from AGN/LINERs by their lower 
[Nil] 6583/Ha and [Om] 5007/H/3 ratios. Among the lines 



proposed to discriminate between the excitation mecha- 
nisms, we plot those from the theoretical modelling of Kew- 
ley et al. (2001) and Stasiriska et al. (2006), and the empiri- 
cal line of Kauffmann et al. (2003), derived from SDSS galax- 
ies. The Shapley galaxies populate both the star-forming and 
the AGN loci. In classifying the emission source, we adopt 
the Stasihska et al. criterion, which is the 'strictest', in the 
sense of assigning fewest galaxies to the normal star-forming 
class. For 23 galaxies with Ha emission, we do not detect 
all of the other three lines with S/N>3 as required for this 
classification method, but eight of these can be confidently 
assigned to the AGN class based on [Nll]/Ha alone. Using 
these criteria we classify 34 galaxies as AGN/LINERs (31 
per cent of all the emission-line galaxies), and 61 objects 
(55 per cent) as pure star-forming galaxies. The remaining 
15 unclassified galaxies have S/N<3 in one or more of the 
lines used, and ambiguous [Nil] /Ha ratio. 

The two categories of emission-line galaxy are dif- 
ferently distributed within the luminosity range spanned 
by our sample. In the upper right panel of Figure [7] we 
show the fraction of galaxies (within the NFPS-selected, 
supercluster-member sample) which have Ha emission, sep- 
arated by BPT-classification, as a function of i?-band lu- 
minosity. Each luminosity bin represents ~70 galaxies. The 
overall emission-line fraction is approximately constant at 
~30 per cent, but the mix between star-forming galaxies 
and AGN/LINERs changes dramatically as a function of 
luminosity. Most of the emission population in the highest- 
luminosity bin are AGN, while star-forming objects predom- 
inate among the faint emission-line galaxies. The break be- 
tween these regimes occurs at Mr ~ —20.5, corresponding 
to luminosity ~ 1O 1O L0 or a stellar mass ~ 10 10 - 4 M Q . This 
transition mass is in excellent agreement with the result ob- 
tained by Kauffmann et al. (2003) for a much larger sample 
of galaxies from SDSS. 

The above analysis used a sample which includes blue 
supercluster members. The lower panels of Figure [7] show 
the same diagrams after imposing an additional selection for 
only the red-sequence galaxies (B — R within 0.15 of the fit in 
Figure [T] 291 galaxies). The overall emission fraction on the 
red sequence is 20 per cent. As expected, the star-forming 



© 2007 RAS, MNRAS 000,111-1201 



12 Russell J. Smith et al. 





Figure 7. Upper left: Emission-line classification diagram, following Baldwin, Phillips & Terlevich (1981), showing the 32 per cent of 
our Shapley sample with measured Ho. The discriminatory lines are from from Kewley et al. (2001) (dotted), Kauffmann et al. (2003) 
(dashed) and Stasihska et al. (2006) (solid). Circular symbols show galaxies from this paper, having sufficient data to classify the likely 
excitation mechanism. Filled circles are galaxies with AGN or LINER emission, according to the Stasinska et al. (2006) criterion, or to 
a simple cut at [Nil] 6583/Ho>0.4. Open circles show normal star-forming galaxies. Crosses indicate galaxies with detected Ha, but 
insufficient S/N in other lines to permit classification. The thick horizontal bar shows the ratio of [Oin] 5007/H/3 = 0.6 -1 often used 
in correcting nebular emission in elliptical galaxy studies. Upper right panel: fraction of AGN/LINERs, star-forming galaxies, and total 
fraction with Ha emission, as a function of total ij-band magnitude. Lower panels: equivalent diagrams for the 20 per cent of Shapley 
red-sequence galaxies with Ha emission. 



branch is much weaker for the red galaxies. The remaining 
galaxies nominally in the star-forming part of the diagram 
appear to form a continuous sequence with the AGN/LINER 
branch (as also seen in the SDSS sample of Yan et al. 2006). 
The BPT-plane sequence is correlated with luminosity: the 
median [Oiii]/H/3 for red-sequence galaxies with Ha emis- 



sion declines by a factor of ten (from 2.1 to 0.24), from 
the most to least luminous galaxy bins (cf. the traditional 
ratio ~ 1.7 = 0.6 -1 ). By contrast the systematic change 
in Ha/H/3 is less than a factor of two. The emission frac- 
tion in our red-sequence sample declines to ~10 per cent for 
fainter galaxies, with the AGN/LINERs cutting off around 
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Figure 8. Distribution of Ha emission equivalent width for the 
sample galaxies, indicating the imposed limit of 0.5 A. Galaxies 
from our primary NFPS-selected sample are shown as the black 
histogram, while the grey section indicates the additional objects 
from 2MASS. The plot limits are such that most of the emission 
line galaxies lie outside of the range shown. 

Mr « —20.5 as before. However, given the systematic trend 
along the sequence with luminosity, it is clear that the AGN 
threshold luminosity depends quite sensitively on the classi- 
fication line adopted (e.g. using the Kewley et al. line would 
result in a cut-off at higher luminosity). Moreover, the ap- 
parently continuous distribution in the BPT plane suggests 
that imposing a classification may even not be justified at 
all, at least for red galaxies. 

To summarize, our analysis reconfirms that the use of 
a constant factor relative to [Oiii] is not a viable means 
of either correcting for or selecting against nebular emis- 
sion contamination. We have shown that most galaxies with 
contaminated H/3 can be readily identified from their Ha 
emission. In the present paper we opt not to make an emis- 
sion correction to the absorption indices, but rather to use 
the Ha measurments to exclude emitting galaxies from the 
sample entirely. 



4 THE FINAL SAMPLE 

This section defines explicitly the sample of low-emission 
supercluster members, and tabulates the data for these ob- 
jects, as analysed in Section [5] 

Recall that after removing galaxies whose entire spec- 
tra were rejected (stars, superpositions, 'fringed' spectra), 
the sample comprises 540 galaxies, with 403 selected from 
the NFPS imaging and 138 from 2MASS. The sample to be 
analysed in this paper will be the subset of NFPS-selected 
objects having very low levels of nebular emission, and hav- 
ing redshifts compatible with membership of the Shapley 
Supercluster. Table Q] shows how the various selection crite- 
ria affect the sample size. 

The emission selection is based on Ha as justified in 
Section [3] In defining the low-emission sample, we use the 
same cut of 0.5 A in EW(Ha), as used in the line-ratio anal- 
ysis. Figure [8] demonstrates that this limit is effective in 
separating the tail of low-emission galaxies from the bulk 
of apparently emission-free objects. Note that the Ha cut 
is equivalent to 0.1-0.2 A in H/?, and is thus more restric- 
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Figure 9. Rcdshift distribution of the observed targets. The black 
histogram shows the NFPS-selected objects, while the grey indi- 
cates the contribution from the 2MASS-selected targets. Galax- 
ies with emission lines are included in the histogram. A total of 
57 galaxies lie beyond z = 0.1 and are not shown. The dashed 
vertical lines indicate the range adopted to define supercluster 
membership. 

tive than the selection at H/3=0.6 A imposed in the NFPS 
analyses (Nelan et al. 2005; Smith et al. 2006). 

The redshift distribution of sample galaxies is shown in 
Figure [9] After an iterative 3cr clipping, the redshift dis- 
tribution has median 14451 km s _1 , with a dispersion of 
927kms~\ A total of 448 observed galaxies (342 NFPS 
plus 106 2MASS) lie within the 3cr range, i.e. 11670- 
17233 kms -1 , which we adopt henceforth as the super- 
cluster redshift limits. In the foreground, there are two 
galaxies with redshifts — 4000kms _1 , characteristic of the 
Hydra-Centaurus Supercluster. A few galaxies have veloci- 
ties — 11000 kms -1 , and are excluded by our limits. (These 
are probably associated with a foreground structure includ- 
ing the rich cluster Abell 3571). A tail of background galax- 
ies extends to z ~ 0.3, with a number of groups appar- 
ent, e.g. at ~23000kms-\ -27000 kms" 1 , -39000 kms -1 , 
-43000 km s" 1 , and -52000 km s" 1 . Several of these fea- 
tures have already been noted from wide-field redshift sur- 
veys of the region (e.g. Proust et al. 2006). 

The final sample analysed in this paper comprises 232 
NFPS-selected galaxies meeting both the emission and the 
redshift criteria, of which 198 have measured velocity disper- 
sions (i.e. are not 'unresolved'). The 34 unresolved galaxies 
would have an average of a = 50 kms -1 , if at the median 
of the Mr — a relationship. However, it is likely that they 
are unresolved precisely because their a is lower than the 
average for their luminosity. 

Referring to the solid points in Figure [T] we note that 
although no explicit colour selection has been imposed, the 
redshift and emission cuts result in a sample which is effec- 
tively limited to the cluster red sequence. A robust fit to 
the B — R colours yields a slope of 0.043 ± 0.005 per mag- 
nitude in R, with a (galaxy-by-galaxy) scatter of 0.07 mag. 
The culled final sample corresponds to —35% of the NFPS 
R < 18 input catalogue. Assuming that the observed and 
unobserved galaxies have similar redshift and emission-line 
distributions, the final sample covers —60% of all passive 
supercluster members within the NFPS imaging region. 

The median measured velocity dispersion is a — 
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104 km s 1 , and the 90 per cent range is 39-228 km s 1 . 
The median error in the HgF index is 0.10 A (0.20 A) for 
velocity dispersions above (below) lOOkrns -1 . Since the 
Thomas et al. (2004) models yield an age-sensitivity of 
d(HgF)/(9(log Age) « —2.6, these values correspond to typ- 
ical formal errors of ~9 per cent (~20 per cent) in age. 

The fully-corrected line-strength data for the final sam- 
ple are presented in Tables [2}j4] 



5 RESULTS 

In this section, we determine the scaling of each of the Lick 
line-strength indices as a function of velocity dispersion, and 
use these to infer average scaling relations for age, metallic- 
ity and a/Fe ratios through comparison to the models of 
Thomas et al. (2003, 2004). The method is essentially that 
of Nelan et al. (2005). 



5.1 Index a relations 

For this analysis, we use the sample of low-emission super- 
cluster members, tabulated in the previous section. The 'un- 
resolved' galaxies (i.e. those which do not have measured a) 
are, of course, not used in the index— a fits, leaving a sam- 
ple of 198 galaxies. Moreover, on an index-by-index basis, 
we remove galaxies with errors greater than four times the 
median (typically only three or four galaxies). The fits were 
weighted according to the error in the line-strengths (ve- 
locity dispersion errors are neglected). An iterative scheme 
was used, in which galaxies deviating from the relation by 
more than four times the intrinsic scatter (see below), and 
by more than twice their measurement error, were rejected 
in each step. However, our results are not very sensitive to 
whether this outlier rejection is applied or not. The total 
and intrinsic scatter around the relation was derived using 
outlier-robust methods: <r to t is the standard deviation of a 
Gaussian distribution having the same interquartile range 
as the residuals from the index— a relation; <Tmt is the ad- 
ditional Gaussian scatter required, in quadrature with the 
measurement errors, to reproduce the observed interquar- 
tile range. For comparison to earlier work on more luminous 
galaxy samples, the relations were also derived for the 104 
galaxies with a > 100 kms" 1 . We have also fit index— a re- 
lations using the traditional VBCs, i.e. a single additive or 
multiplicative correction irrespective applied to all galaxies 
of given a. The new VBCs affect the slopes by <3%, with 
the exceptions of: HdA 7% shallower with new correction, 
Ca4227 11% shallower and G4300 8% steeper. 

The index— cr relations are presented in Figure [lOl and 
summarized in Table [S] Comparison of our a > lOOkrns -1 
slopes with the relations tabulated by Nelan et al. reveals 
that the slopes are typically ~20 per cent flatter in the 
present work. The only indices not consistent with this factor 
are CN1 and Ca4227. The difference in slope is in the sense 
expected from the change in velocity dispersion scale, but 
larger than suggested by the a comparison in Section 12.31 
For Ca4227, the resolution corrections applied by Nelan et 
al. were erroneous, generating a sharp increase in the index 
at high a and consequently a much steeper slope. Ca4227 
is the only case where this effect is non-negligible, since its 
pass-bands are much narrower than any other of the other 



indices. (Note that Ca4227 was not used in the Nelan et 
al. stellar population analysis.) For CN1, where there is ap- 
parently some curvature in the relation, our slope for the 
a > 100 km s -1 set is slightly steeper than in Nelan et al., 
while our slope for the full sample is flatter. 



5.2 The stellar population scaling relations 

Consideration of the index— cr slopes, and comparison to 
stellar population models, provides a means to determine the 
average scaling of age and metal abundances as a function of 
mass (e.g. Nelan et al. 2005; Clemens et al. 2006). Although 
using only the slopes forces the scaling relations to a very 
simple (and perhaps inaccurate) form, this is an explicitly 
relative method, and therefore is insensitive to zero-point 
offsets in either the data or the models. Further discussion 
of the merits and limitations of the method can be found in 
Nelan et al. and in Smith et al. (2006). 

From the Thomas et al. (2003, 2004) model grids, we de- 
termine the responses of each index to the three population 
parameters by fitting a linear relation I = J?A go log Age + 
-Rz/h[Z/H] + _R a / Fc [a/Fe], for models with age 3-15 Gyr, 
[a/Fe]=0.0-0.5 and [Z/H]=-0.33, 0.00, +0.35, +0.67. This is 
a slight improvement on the original method where fiducial 
values were assumed for two parameters when deriving the 
index response of the third parameter. For reference, we pro- 
vide the new responses in Table [S] 

To recover the scaling relations of the three stellar pop- 
ulation parameters, we consider the observed slopes for nine 
of our index-a relations (HdF, HgF, Hbeta, CN1, Fe4383, 
Mgb5177, Fe4668, Fe5015 and Fe5406). These are chosen to 
be mainly non-redundant (e.g. only one H7 index is used), 
to avoid indices which are contaminated by the 5577 A sky 
line in many cases (i.e. Fe5270 and Fe5335), and to avoid in- 
dices with strong dependence on Ca abundancfl The G4300 
index is not used in this analysis since it disagrees strongly 
with the other indices, and severely degrades the fit qual- 
ity if included. The G4300-a relation is much flatter than 
expected from the other indices, perhaps due to the high- 
metallicity saturation mentioned by Schiavon (2007). 

The scaling relations derived through the comparison of 
measured slopes with the model-derived responses are: 



A 

Age oc a 



\ Z/H 



*/Fe 



These scalings can be used to predict the index— a slopes 
for a wider range of indices, as given in Table [6] which can 
be directly compared to the observed slopes. The quoted er- 
rors are the formal uncertainty in the fit, and do not include 
modelling uncertainties such as the scatter when different 
indices are used for the fit. The latter can be estimated con- 
servatively by repeating the exercise for all non-degenerate 
"triplets" chosen from the above indices. We impose the con- 
dition that each triplet must include one Balmer line, one 
index from the set (Mgb5177, CN1, Fe4668) where R a/Fc 
has the same sign as Rz/n and one from the set (Fe4383, 
Fe5015, Fe5406) where [a/Fe] and [Z/H] responses differ in 
sign. The scatter in scaling relation exponents, computed 

2 Ca is nominally an a element, and enhanced in the default 
Thomas et al. (2003, 2004) models, but is observed to track Fe 
rather than Mg (Cenarro et al. 2004, and references therein). 
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Table 1. Sample size after application of various selection criteria. 



Selection criteria (cumulative) 


Cut applied 




Number of galaxies 








NFPS-selccted 2MASS-selected 


Total 


all observed 




416 


149 


565 


uncorrupted galaxy spectra 


See Section l2~2l 


402 


138 


540 


A range for a 


cz < 22000 km s _1 


350 


111 


461 


low nebular emission 


EW(Hcv)<0.5A 


232 


50 


282 


Shapley member 


11670 kms-k cz <17233kms- 1 


232 


48 


280 


resolved cr 


See Section [231 


198 


47 


245 



Table 2. Basic data and principal line-strength indices for the sample galaxies. The table includes data for the 232 galaxies in the 
NFPS-selected, supercluster-member and emission-free sample as analysed in this paper. The galaxy positions can be obtained from 
their identification numbers. S/N is the signal-to-noise ratio measured over a rest frame interval of 4500— 5500 A. The redshift cz is in 
the heliocentric frame. Columns R and B are total magnitudes from the NFPS photometry. Velocity dispersions are given as log a (with 
a in kms -1 ), when measured. Galaxies with missing data in this column are those for which the dispersion was 'unresolved'. This table 
includes three non-degenerate absorption line indices Mgb5177, Fe5015 and HgF which can be used to invert stellar population models 
(see Paper II). Indices have been corrected to the Lick resolution and to zero intrinsic velocity broadening. Velocity dispersions and line 
indices are as observed through a 2arcsec diameter aperture, corresponding to 1.9 kpc in our adopted cosmology. Additional indices are 
provided in Tables l3l and l4l 



Galaxy ID 




S/N 


CZ 


R 


B 


logo- 


Mgb5177 


Fe5015 


HgF 


NFPJ132328.9- 


-314242 


60 


13940 


16.127 


17.964 


2.067 ±0.019 


4.107 ±0.124 


4.886 ± 0.210 


-1.384 ±0.144 


NFPJ132330.8- 


-314935 


93 


15061 


15.467 


17.185 


2.014 ±0.012 


3.494 ± 0.074 


4.726 ± 0.137 


-0.656 ± 0.089 


NFPJ132335.5- 


-315201 


40 


14870 


16.739 


18.577 


2.005 ±0.029 


3.410 ±0.180 


5.576 ± 0.307 


-1.277 ±0.206 


NFPJ132337.1- 


-315047 


44 


14337 


16.163 


17.984 


2.129 ±0.017 


3.812 ±0.152 


4.021 ± 0.273 


-1.144 ±0.170 


NFPJ132345.0- 


-314230 


54 


14999 


16.292 


18.118 


1.843 ±0.039 


2.854 ±0.127 


4.332 ± 0.229 


-0.031 ± 0.152 


NFPJ132348.3- 


-314953 


38 


15123 


17.240 


19.061 


1.810 ±0.038 


3.649 ±0.179 


4.617 ±0.316 


-0.913 ± 0.240 


NFPJ132350.3- 


-313519 


37 


14034 


17.330 


19.023 




2.611 ±0.196 


4.348 ± 0.317 


0.331 ± 0.235 


NFPJ132355.5- 


-313847 


34 


15248 


17.838 


19.645 


1.861 ±0.046 


3.780 ±0.193 


4.969 ± 0.344 


-1.721 ± 0.301 


NFPJ132406.9- 


-314449 


58 


13600 


16.404 


18.260 


2.041 ±0.019 


3.648 ±0.137 


4.426 ± 0.220 


-1.097 ±0.154 






The 


full content 


of this data table is 


available in the electronic version of the 


paper. 





Table 3. Supplementary line-strength index data for the sample galaxies. 



Galaxy ID HdA HdF CN1 CN2 Ca4227 G4300 



NFPJ132328.9- 


-314242 


-1.43 ±0.27 


0.56 ±0.19 


0.031 ± 0.007 


0. 


.065 ± 0.009 


1.32 ± 0.10 


5. 


.60 ±0.22 


NFPJ132330.8- 


-314935 


-0.77 ±0.15 


1.09 ±0.10 


0.016 ±0.004 





.045 ± 0.005 


1.06 ± 0.06 


4. 


.62 ±0.13 


NFPJ132335.5- 


-315201 


-1.40 ±0.38 


0.49 ±0.26 


0.038 ± 0.010 





.077 ±0.012 


1.28 ± 0.14 


5. 


.39 ±0.30 


NFPJ132337.1- 


-315047 


-1.20 ±0.30 


0.35 ±0.21 


0.030 ± 0.008 





.060 ± 0.010 


0.96 ± 0.12 


4. 


.94 ±0.26 


NFPJ132345.0- 


-314230 


0.68 ± 0.24 


1.12 ±0.18 


-0.034 ± 0.007 - 


-0 


.010 ± 0.008 


1.07 ±0.09 


4. 


.50 ±0.23 


NFPJ132348.3- 


-314953 


-1.95 ± 0.43 


0.32 ±0.30 


0.018 ±0.011 


0. 


.056 ± 0.014 


1.33 ± 0.15 


5. 


.89 ±0.36 


NFPJ132355.5- 


-313847 


-0.64 ±0.55 


0.88 ±0.39 


-0.020 ±0.015 - 


-0 


.000 ± 0.018 


1.50 ±0.19 


5. 


.36 ±0.44 


NFPJ132406.9- 


-314449 


-1.12 ± 0.28 


0.36 ±0.20 


0.016 ±0.008 





.041 ± 0.009 


1.24 ±0.11 


1 


.92 ±0.24 


NFPJ132412. 7-314658 


-0.05 ± 0.33 


1.04 ±0.24 


0.006 ± 0.009 





.037 ±0.011 


1.12 ± 0.12 


4. 


.20 ±0.29 




The 


full content of th 


is data table is 


available in the electron 


ic 


version of the 


paper. 







Table 4. Supplementary line-strength index data for the sample galaxies. 



Galaxy ID 




HgA 


Fe4383 


Ca4455 




Fe4531 


Fe4668 


Hbcta 




Fe5406 


NFPJ132328.9- 


-314242 


-5.68 ±0.26 


4.67 ±0.28 


1.35 ±0.11 


3. 


.57 ±0.22 


5.91 ±0.30 


2.06 ± 0. 


11 


1 


.91 ±0.11 


NFPJ132330.8- 


-314935 


-4.30 ±0.16 


4.58 ±0.17 


1.17 ±0.07 


3 


.34 ±0.16 


5.55 ± 0.19 


2.02 ±0. 


07 


1 


.46 ±0.08 


NFPJ132335.5- 


-315201 


-5.84 ±0.38 


5.31 ±0.40 


1.41 ±0.15 


3. 


.20 ±0.34 


5.33 ± 0.44 


1.80 ±0. 


16 


1 


.82 ±0.18 


NFPJ132337.1- 


-315047 


-5.43 ±0.31 


4.45 ±0.35 


0.95 ±0.14 


2 


.92 ±0.27 


4.60 ±0.38 


1.60 ±0. 


14 


1 


.55 ±0.15 


NFPJ132345.0- 


-314230 


-3.18 ±0.28 


3.89 ±0.31 


0.97 ±0.11 


2 


.84 ±0.27 


4.29 ± 0.35 


2.25 ±0. 


12 


1 


.21 ±0.13 


NFPJ132348.3- 


-314953 


-5.28 ±0.44 


5.18 ±0.44 


1.22 ±0.16 


2 


.70 ±0.39 


4.72 ± 0.48 


1.92 ±0. 


17 


1 


.65 ±0.18 


NFPJ132355.5- 


-313847 


-5.78 ±0.54 


4.61 ±0.52 


1.75 ±0.19 


3 


.27 ±0.43 


5.14 ±0.52 


1.86 ±0. 


19 


1 


.57 ±0.19 


NFPJ132406.9- 


-314449 


-4.77 ±0.27 


3.86 ±0.32 


1.02 ±0.12 


2 


.79 ±0.23 


4.37 ±0.31 


1.75 ±0. 


12 


1 


.56 ±0.11 


NFPJ132412.7- 


-314658 


-4.66 ±0.33 


4.84 ±0.36 


1.24 ±0.13 


3. 


.23 ±0.28 


5.66 ± 0.40 


2.16 ±0. 


14 


1 


.56 ±0.15 



The full content of this data table is available in the electronic version of the paper. 
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log(a) log(a) log(a) log(a) 

Figure 10. Indcx-a relations. The lines indicate the fits to all of the data (solid) and to the a > 100 km s" 1 subset (dashed). The objects 
with unresolved velocity dispersions (to the left of the dotted line) are not used when fitting the relations. Open symbols represent these 
and also other galaxies excluded from the fit by an iterative outlier rejection scheme. 



over all 27 of these triplets is 0.10 for age, 0.07 for [Z/H] 
and 0.06 for [a/Fe], which may be taken as representative 
of the systematic errors. As a further test, we can use all six 
metal lines plus each Balmer line in turn. In this case, the 
age exponents are 0.59 using HdF, 0.44 using HgF and 0.54 
using Hbeta, with random errors 0.08 in each case. Fitting 
the scaling relations to the index— a slopes derived using the 
"traditional" VBCs (i.e. a common additive or multiplica- 
tive correction function for all spectral types), the exponents 
change by less than 0.01 relative to our default case. Thus we 
do not support the claim by Kelson et al. (2006, their Figure 
10) that subtle differences in VBC method can dramatically 
affect the derived relations. 

Applying the same method to our index— a slopes de- 
rived from the high-mass (a > 100 km s -1 ) subset, we ob- 
tain: 



» 0.64±0.12 r, m 0.38±0.09 /-ri 0.36±0.07 

Age tx a , Z/H oc a , a/be oc a 



The high-cr age and [Z/H] trends are consistent with the full 
sample within la, while there is a ~2<r change in slope for 
[a/Fe]. The latter suggests a flattening of the [a/Fe]— a re- 
lation towards low masses. The difference is driven by the 
a-sensitive Mgb5177 and CN1 indices. In particular a curva- 
ture, or change in slope around log a = 2, is visually appre- 
ciable in the CN index panels of Figure fTTJl Similar behaviour 
can be seen in the data of Sanchez-Blazquez et al. (2006a, 
figure 6). 

The recovered parameter scaling relations from the 
above analysis are summarized in the upper section of Ta- 
ble 
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Table 5. Summary of the index— a relations of Figure lTOl For each index, the table gives the zero-point and slope of the fitted relationship, 
together with robust estimates of the total and intrinsic scatter, trtot ar >d fint- Results are given both for the full sample and for the 
subset with a > 100 km s . 



Index All galaxies (N=198) a > 100 km s" 1 (N=104) 

zero-point slope trtot o"mt zero-point slope otot 0i„t 
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Table 6. Linear responses for indices derived from model grids. 

ar is the change in the index for a decade change in par, where 
par is Age, Z/H or a/Fe. The responses are derived from multi- 
variate fits to the Thomas et al. (2003, 2004) models with age 
3-15 Gyr, [a/Fe]=0.0-0.5, and [Z/H]=-0.33, 0.00, +0.35, +0.67. 
The final column uses the scaling relations obtained in Section [5] 
to predict the slope of the index— a slope. 



index 


-^Agc 


^?Z/H 


R-a/Fe 


predicted slope 


HdA 


-3.816 


-3.614 


+4.707 


-2.130 


HdF 


-1.807 


-1.278 


+1.954 


-0.925 


CN1 


+0.131 


+0.199 


+0.059 


+0.149 


CN2 


+0.133 


+0.208 


+0.067 


+0.155 


Ca4227 


+0.877 


+1.193 


+0.258 


+0.921 


G4300 


+1.513 


+1.876 


+1.482 


+1.765 


HgA 


-4.383 


-3.974 


+3.414 


-2.845 


HgF 


-2.630 


-2.116 


+0.999 


-1.857 


Fe4383 


+1.976 


+3.753 


-4.099 


+1.361 


Ca4455 


+0.582 


+1.224 


+0.062 


+0.733 


Fe4531 


+0.865 


+1.559 


-0.941 


+0.764 


Fe4668 


+1.884 


+7.929 


+0.503 


+3.791 


Hbeta 


-1.153 


-0.585 


+0.285 


-0.733 


Fe5015 


+0.826 


+2.713 


-1.422 


+1.025 


Mgb5177 


+1.576 


+2.971 


+2.076 


+2.307 


Fe5406 


+0.448 


+1.158 


-1.182 


+0.355 



6 DISCUSSION 

In the lower part of Table[7] we tabulate for comparison some 
results from other recent studies. The analyses were based 
either on the index— a slopes method, or on explicit inversion 
of model grids. Most of the studies used the Thomas et al. 
(2003, 2004) models. Exceptions are Graves et al. (2007) 
(based on models of Schiavon 2007), Sanchez-Blazquez et 
al. (2006b) (based on Vazdekis et al. 2007) and Clemens et 
al. 2006 (based on Annibali et al. 2005). Before discussing 
the comparisons we again emphasize that our sample spans 
a factor of ~6 in velocity dispersion, compared to a factor 
^2-3 in the other work. 



The most comparable studies, using similar galaxy se- 
lection criteria, and similar methods, are those from the 
NFPS (Nelan et al. 2005; Smith et al. 2006). The two NFPS 
analyses were based on slightly different samples and data 
treatment, and the differences between their results proba- 
bly reflect the level of systematic errors in the method. Our 
results from this paper, for the high-rr subset, are within 
the range spanned by the two NFPS analyses. For the full 
sample, i.e. extending into the low-cr regime, the trends are 
slightly natter than in NFPS, especially for [a/Fe]. A possi- 
ble change in some of the index— a slopes has already been 
commented on, but the change in scaling relations may also 
reflect the difference in a scales discussed in Section \2. 31 Fi- 
nally, the inclusion of the low mass galaxies may be compro- 
mising a key assumption of the slopes method, namely that 
the grids are linear and parallel over the index range spanned 
by the data. In Paper II, we will bypass this requirement by 
explicitly inverting the model grids on a galaxy-by-galaxy 
basis. 

Among other work, results based on samples from SDSS 
are similar to NFPS and to the present survey, in targeting 
red-sequence or early-type galaxies without detailed visual 
inspection. The work of Bernardi et al. (2006) agrees fairly 
well with our results and with NFPS, although all three pa- 
rameter slopes are rather steeper (we have quoted their re- 
sults using HgF as the age indicator; their age slope is much 
steeper when Hbeta is used.) The work of Clemens et al. 
(2006) agrees in age slope, but yields much steeper [Z/H] and 
[«/Fe] scalings, which may result from the different models 
they use, tracking the C abundance separately from Mg and 
Fe. In comparing between NFPS, SDSS and our new sam- 
ple, note that the mass baseline covered by SDSS is a factor 
of ~3 narrower. Most recently, Graves et al. (2007) have ob- 
tained a relatively shallow age slope and steeper metallicity 
slope from a colour- and emission-selected SDSS sample. 

The compilation of Thomas et al. (2005) yields scal- 
ings close to those of NFPS, when a fit is made to their 
whole sample, rather than excluding young low-mass galax- 
ies in advance as in their paper. The work of Caldwell et al. 
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(2003), which was among the first to obtain age estimates in 
the a < 100 km s -1 regime, suggests a steeper age slope than 
recovered here, and a comparable metallicity trend. (Cald- 
well et al. correct for the a-enhancement effect in deriving 
the metallicities, but do not report [a/Fe] measurements for 
their sample.) The recent study by Sanchez-Blazquez et al. 
(2006b), on the other hand, yields a very flat age— a relation 
in the Coma cluster (from ~20 galaxies), although a fairly 
steep relation appears to hold in their Virgo and field/group 
galaxies. The discrepancy between y Sanchez-Blazquez et al. 
and the other results is striking, and remains to be under- 
stood. In particular, results from larger samples in Coma 
(e.g. Matkovic & Guzman 2005) should help to resolve this 
issue. 

On balance, most of the recent line-strength-based stud- 
ies concur on a significant age trend with cr, such that while 
the stars in the most massive galaxies are almost as old as 
the universe, the low-cr galaxies on average continued to form 
stars into the z < 1 era. For example, taking Agecx a 0,5 , 
and forcing the age at 300kms _1 to be 12Gyr, our trend 
requires an average age of ~7Gyr at a — 100 km s^ 1 , and 
~5Gyr at a = 50kms _1 . The corresponding redshift of the 
last star-formation epoch is z « 0.8 for a — 100 km s , and 
z ~ 0.5 for a = 50kms _1 . The implications of these results 
for high-redshift observations have been discussed elsewhere 
(e.g. Smith 2005; Nelan et al. 2005), but can be updated 
here. Specifically, some 50 per cent of the faint galaxies, now 
residing on the red sequence ~2 magnitudes below L* , were 
star-forming galaxies at z ~ 0.5. The observed age-mass re- 
lation therefore suggests that a deficit of faint red-sequence 
galaxies should be observed in high redshift clusters. Early 
claims for such a deficit (De Lucia et al. 2004; Kodama et 
al. 2004) have been supported by more extensive observa- 
tions (De Lucia et al. 2007; Muzzin et al. 2007; Stott et al. 
2007), although contrary results have been obtained (An- 
dreon 2006). De Lucia et al. define a giant-to-dwarf ratio, 
G/D, on the red sequence of 20 clusters, based on number 
of red galaxies in luminosity intervals corrected for passive 
evolution. They report an increase of G/D by a factor of ~3 
between z — and z — 0.75. We have conducted some simple 
Monte-Carlo experiments using the red sequence luminosity 
function and the luminosity-a correlation, to translate our 
age-cr relation into estimates for the evolution in G/D. These 
tests yield G/D increases by a factor of 2-5 at z — 0.75, for 
Agecx cr ' 5 , where the range includes different choices of lu- 
minosity function parameters, different normalisations of the 
age-cr relation etc. We conclude that the ages measured in 
present day red-sequence galaxies not only require a deficit 
of faint red galaxies in high redshift clusters as observed, but 
also agree quantitatively in the level of of G/D evolution out 
to z = 0.75. 

As a closing caveat, we note that although the age-cr 
relation has become a standard diagnostic for downsizing, it 
is usually applied to magnitude-limited samples. At a given 
cr, galaxies which are younger than average for their mass 
will be brighter than average, and so preferentially included 
in such samples, with the result of steepening the recovered 
trends. The impact of the bias depends on the limiting lu- 
minosity reached by the sample, and is particularly severe 
for surveys such as SDSS which cover only the exponential 
cut-off in the cr distribution. This effect was discussed by 
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Figure 11. Consistency between present day red-sequence ages 
and the observed redshift evolution of the giant-to-dwarf galaxy 
ratio. The solid line shows the expected evolution of G/D (defined 
as in De Lucia et al. 2007), if the red-sequence galaxies have Age oc 
cr 0,5 as found in this paper, and if the most massive ellipticals 
formed their stars at z ~ 3. The calculation assumes generic 
parameters for the luminosity function and the Mr — a relation. 
Dotted lines repeat the calculation for exponents 0.7, 0.6, 0.4, 0.3 
(left to right). Points with error bars represent measurements of 
G/D in Coma and in the ESO Distant Cluster Survey (De Lucia 
et al. 2007). 

Kelson et al. (2006), and is investigated in greater detail in 
Paper II. 



7 CONCLUSIONS 

This paper has presented the first phase of a survey of galaxy 
properties in the Shapley Supercluster. We have described 
the processing of spectra obtained from ~8 hr integrations 
with AAOmega. As well as redshifts and velocity disper- 
sions, we have measured emission line equivalent widths and 
absorption line indices. The latter incorporate an improved 
treatment of corrections for instrumental resolution and ve- 
locity broadening. 

We find a 32 per cent incidence of Ha emission in the 
supercluster sample. For red-sequence galaxies, emission is 
detected in 20 per cent, with line ratios varying systemati- 
cally from AGN-like at high luminosity to being more consis- 
tent with normal star formation at low luminosity. Defining 
a low-emission subsample based on the Ha measurements, 
we have presented the index— a relations for a range of Lick 
line-strength indices. We show that the index— a relations 
and the physical scaling relations continue into the low-mass 
regime with little change in slope in most cases. Considering 
the slopes of these relations, in the context of the Thomas 
et al. (2003, 2004) stellar population models, we derived 
estimates for the scaling relations of age, total metallicity 
and a-element abundance ratio. All three parameters are 
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Table 7. Comparison of parameter-vs-cr scaling relations from different index sets in this work, and with recent results from large samples 
(^Vgal ^ 100). To emphasize the different mass ranges covered by different studies, we indicate the median velocity dispersion, {a), and 
the 90 per cent range (estimated in some cases). The column headed 'span' is the ratio of a spanned by the 90 per cent range. 



Reference 


source 


Method 


90% range cr 


<^> 


span 


age slope 


[Z/H] slope 


[cv/Fe] slope 


note 


A lllo W Wl xv 






39- 


-228 


104 


5.8 


n ko + fl f)6 


L/.tJrt _i_ yj.\j'-t 




I 
















± 0.10 


± 0.07 


± 0.06 


2 


This work (cr > 100 km s 




slopes 


104- 


-242 


152 


2.3 


0.64 ±0.12 


0.38 ± 0.09 


0.36 ±0.07 


1 
















±0.23 


±0.15 


±0.13 


2 


Nelan et al. (2005) 


NFPS 


slopes 


74- 


-260 


147 


3.5 


0.59 ±0.13 


0.53 ±0.08 


0.31 ±0.06 


3 


Smith et al. (2006) 


NFPS 


slopes 


74- 


-260 


147 


3.5 


0.72 ±0.14 


0.37 ±0.08 


0.35 ±0.07 


4 


Graves et al. (2007) 


SDSS 


inversion 


100- 


-250 


175 


2.5 


0.35 ±0.03 


0.79 ±0.05 


0.36 ±0.04 


5 


Bernardi et al. (2006) 


SDSS 


inversion 


120- 


-265 


188 


2.2 


0.81 


0.58 


0.39 


6 


Clemens et al. (2006) 


SDSS 


hybrid 


120- 


-270 


180 


2.2 


~0.8 


0.76 


0.74 


7 


Thomas et al. (2005) 


Compilation 


inversion 


106- 


-310 


208 


2.9 


0.48 ±0.21 


0.35 ±0.09 


0.40 ±0.05 


8 


Caldwell et al. (2003) 




inversion 


58- 


270 


153 


4.6 


0.80 


0.32 




9 



Notes: errors are formal uncertainty from fit to multiple indices; systematic errors from dispersion among non-degenerate index triplets (see text); 
^ error includes systematic component; 4 formal error only; ^ Hbeta ages for quiescent sample, formal error only; ^ results using HgF (their Table 9); 
7 metallicity trends reported by Clemens et al., our estimate of age trend from their Figure 10; ^ our fit (error- weighted) to all data from Thomas et al., 
formal error only; metallicity trend as reported, our estimate of age trend from their Table 5. 



found to increase significantly with a, as found also in sam- 
ples limited to high-mass red-sequence galaxies. There is a 
slight tendency towards natter scalings at low a, especially 
in the case of [a/Fe]. Focusing on the age scaling, our esti- 
mated relation Age oc a ' 5 is quantitatively consistent with 
measurements of dwarf-deficient red sequences at redshifts 
z = 0.4 - 0.8 e.g. by De Lucia et al. (2007). 

Our study is unique in including a large number of low- 
luminosity galaxies, with sufficient signal-to-noise to obtain 
meaningful age and abundance estimates on a galaxy-by- 
galaxy basis. This analysis is presented in Paper II (Smith et 
al. 2007). Future papers will combine the AAOmega spectro- 
scopic data with imaging from the Shapley Optical Survey 
(Mercurio et al. 2006; Haines et al. 2006) to address stel- 
lar population variations as a function of morphology and 
environment. Observational extensions to the spectroscopic 
survey are also underway. 
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